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Abstract. The algebraic structure and the spectral properties of a special class of multi- 
£Nj . component NLS equations, related to the symmetric spaces of BD.I-type are analyzed. 

The focus of the study is on the spectral theory of the relevant Lax operators for different 
fundamental representations of the underlying simple Lie algebra q. Special attention is paid 
^0 ' to the structure of the dressing factors in spinor representation of the orthogonal simple Lie 

£H , algebras of B r ~ so(2r + 1, C) type. 
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> • 

o 

^ ! 1 Introduction 

. The nonlinear Schrodinger equation [49|, [2] 

o ' 

iqt + Qxx + 2|g| 2 <? = 0, q = q(x,t) 

has natural multi-component generalizations. The first multi-component NLS type model with 
applications to physics is the so-called vector NLS equation (Manakov model) |40[ [2] : 



iv t + v xx + 2(v T , v)v = 0, v 



/ vi{x,t) 



\v n (x,t) 



Here v is an n-component complex- valued vector and (•, •) is the standard scalar product. All 
these models appeared to be integrable by the inverse scattering method [H \2\ Wf \ \TT \ [6| IT] |31~ | [36] . 

Later on, the applications of the differential geometric and Lie algebraic methods to soliton 
type equations [101 SU LH [211 [23l [39l [Ml E21 [26l [371 ESI O E] (for a detailed review see e.g. |M]) 
has lead to the discovery of a close relationship between the multi-component (matrix) NLS 
equations and the homogeneous and symmetric spaces [12]. It was shown that the integrable 
MNLS type models have Lax representation with the generalized Zakharov-Shabat system as 
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the Lax operator: 

Lip(x,t,X) =i^ + (Q(x,t) - XJ)tp(x,t,X) = 0, (1.1) 
ax 

where J is a constant element of the Cartan subalgebra f) C of the simple Lie algebra g and 
Q(x,t) = [J,Q{x,t)\ G g/f). In other words Q(x,t) belongs to the co-adjoint orbit M j of g 
passing through J. The Hermitian symmetric spaces, compatible with the NLS dispersion law, 
are labelled in Cartan classification by A. Ill, C.I, D.III and BD.I, see [5*1 133] . 

In what follows we will assume that the reader is familiar with the theory of simple Lie 
algebras and their representations. The choice of J determines the dimension of Mj which can 
be viewed as the phase space of the relevant nonlinear evolution equations (NLEE). It is equal 
to the number of roots of g such that a(J) ^ 0. Taking into account that if a is a root, then 
—a is also a root of g, so dimMj is always even. 

The interpretation of the ISM as a generalized Fourier transforms and the expansion over 
the so-called 'squared' solutions (see |14| for regular and [20} [22] for non-regular J) are based 
on the spectral theory for Lax operators in the form (jl.ip . This allow one to study all the 
fundamental properties of the corresponding NLEE's: i) the description of the class of NLEE 
related to a given Lax operator L(X) and solvable by the ISM; ii) derivation of the infinite family 
of integrals of motion; and iii) their hierarchy of Hamiltonian structures. 

Recently, it has been shown, that some of these MNLS models describe the dynamics of spinor 
Bose-Einstein condensates in one-dimensional approximation |34j [44] . It also allows an exact 
description of the dynamics and interaction of bright solitons with spin degrees of freedom |38] . 
Matter- wave solitons are expected to be useful in atom laser, atom interferometry and coherent 
atom transport (see e.g. [45] and the references therein). Furthermore, a geometric interpretation 
of the MNLS models describing spinor Bose-Einstein condensates are given in [12]; Darboux 
transformation for this special integrable case is developed in [38]. 

Along with multi-component NLS-type of systems, generalizations for other hierarchies has 
also attracted the interest of the scientific community: The scalar |46] and multi-component 
modified Korteweg-de Vries hierarchies over symmetric spaces [3] have been further studied 
in [121 El [28]. 

It is well known that the Lax representation of the MNLS equations takes the form: 
[L, M] = 

with conveniently chosen operator M, see equation (|2.2|) below. In other words the Lax repre- 
sentation is of pure Lie algebraic form and therefore the form of the MNLS is independent on the 
choice of the representation of the relevant Lie algebra g. That is why in applying the inverse 
scattering method (ISM) until now in solving the direct and the inverse scattering problems 
for L (jl.ip only the typical (lowest dimensional exact) representation of g was used. 

Our aim in this paper is to explore the spectral theory of the Lax operator L (jl.ip for 
different fundamental representations of the underlying simple Lie algebra g. We will see that 
the construction of such important for the scattering theory objects like the fundamental an- 
alytic solutions (FAS) depend crucially on the choice of the representation. This reflects on 
the formulation of the corresponding Riemann-Hilbert problem (RHP) and especially on the 
structure of the so-called dressing factors which allow one to construct the soliton solutions of 
the MNLS equations. In turn the dressing factors determine the structure of the singularities 
of the resolvent of L. In other words one finds the multiplicities of the discrete eigenvalues of L 
and the structure of the corresponding eigensubspaces. 

We will pay special attention to the adjoint representation, which gives the expansion over 
the so-called 'squared' solutions. For MNLS related to the orthogonal simple Lie algebras of B r 
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and D r type, we will outline the spectral properties of the Lax operator in the spinor represen- 
tations. Another important tool is the construction of the minimal sets of scattering data T,, 

1 = 1,2, each of which determines uniquely both the scattering matrix T(A) and the poten- 
tial Q(x). Our remark is that the definition of Tj used in [141 [20} [25| [26] is invariant with 
respect to the choice of the representation. 

The paper is organized as follows: In Section [2j we give some preliminaries about MNLS 
type equations over symmetric spaces of BD.I-type. Here we summarize the well known facts 
about their Lax representations, Jost solutions and the scattering matrix T(A) for the typi- 
cal representation, see j!4[ I16j. Next we outline the construction of the FAS and the relevant 
RHP which they satisfy. They are constructed by using the Gauss decomposition factors of 
the scattering matrix T(A). All our constructions are applied for the class of potentials Q(x) 
that vanish fast enough for x — > dboo. We finish this section by brief exposition of the simplest 
type of dressing factors. We also introduce the minimal sets of scattering data as the minimal 
sets of coefficients Tj, i = 1,2, which determine the Gauss factors of T(A). These coefficient 
are representation independent and therefore Tj determine T(A) and the corresponding poten- 
tial Q{x) in any representation of g. In Section [3] we describe the spectral properties of the Lax 
operator in the typical representation for BD.I symmetric spaces and the effect of the dressing 
on the scattering data. We introduce also the kernel of the resolvent R ± (x, y, A) of L and use 
it to derive the completeness relation for the FAS in the typical representation. This relation 
may also be understood as the spectral decomposition of L. In Section H] we describe the spec- 
tral decomposition of the Lax operator in the adjoint representation: the expansions over the 
'squared' solutions and the generating (recursion) operator. Most of the results here have also 
been known for some time [TJ [T4"l 131] . In particular the coefficients of Tj appear as expansion 
coefficients of the potential Q{x) over the 'squared solutions'. This important fact allows one 
to treat the ISM as a generalized Fourier transform. In the next Section [5] we study the spectral 
properties of the same Lax operator in the spinor representation: starting with the algebraic 
structure of the Gauss factors for the scattering matrix, associated to L(X), the dressing factors, 
etc. Note that the FAS for BD.I-type symmetric spaces are much easier to construct in the 
spinor representation. One can view L sp as Lax operator related to the algebra so(2 r ) with 
additional deep reduction that picks up the spinor representation of g ~ so(2r + 1). In the last 
Section [6] we extend our results to a non- fundamental representation. In order to avoid unneces- 
sary complications we do it on the example of the 112-dimensional representation of so{7) with 
highest weight uj = 3u>3 = |(ei + e2 + 63). The matrix realizations of both adjoint and spinor 
representations for B2 ^ so(5,C) and B7 ~ so(7, C) are presented in Appendices [A) IB1 and ICl 

2 Preliminaries 

We start with some basic facts about the MNLS type models, related to BD.I symmetric 
spaces. Here we will present all result in the typical representation of the corresponding Lie 
algebra g ~ B r ,B r . These models admit a Lax representation in the form: 



Here, in general, J is an element of the corresponding Cartan subalgebra f) and Q(x,t) is an 
off-diagonal matrix, taking values in a simple Lie algebra g. In what follows we assume that 
Q{x) G Mj is a smooth potential vanishing fast enough for x — > ±00. 

Before proceeding further on, we will fix here the notations and the normalization conditions 
for the Cartan-Weyl basis {hk,E a } of g (r = rankg) with a root system A. We introduce 



Lip(x, t, A) = id x ip + (Q(x, t) - XJ)ip(x, t, A) = 0, 

M^(x,t,X) = id t ^+ (V Q (x,t) + XVi(x,t) - X 2 J)^{x,t,X) = 0, 

V!(x,t) = Q(x,t), V (x,t) =iadj 1 ^ + ^[adj 1 Q,Q(x,t)] . 



(2.1) 
(2.2) 
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hk G f), k = 1, . . . , r as the Cartan elements dual to the orthonormal basis {e k } m the root 
space W and the Weyl generators E a , a G A. Their commutation relations are given by [51 133]: 

2 - 

[h k ,E a ] = (a, ejfc)£7 a , [E a , E- a ] = - r }(&, e k )h k , 

(a, a) ^ 



[i£ a , 



N a ,pE a+ p for a + /3 G A, 

for a + /3 0Au{O}. 



Here a = X)fc=i a fc e fc ^ s a ^-dimensional vector dual to J G f) and (•, •) is the scalar product in E r . 
The normalization of the basis is determined by: 

E- a = E , (E_ a ,E a ) = - -, N- a p = —N a p, 

[a, a) 

where N a p = ±(p + 1) and the integer p > is such that a + s(3 G A for all s = 1, . . . ,p, 
a + (p + A and (■, •) is the Killing form of g, see [5j [33]. The root system A of g is 
invariant with respect to the group W Q of Weyl reflections S a , 

Jay = V ~ ~, r«, a £ A. 

[a, a) 

As it was already mentioned in the Introduction the MNLS equations correspond to Lax opera- 
tor ([l.ip with non-regular (constant) Cartan elements J G f). If J is a regular element of the 
Cartan subalgebra of g then ad j has as many different eigenvalues as is the number of the roots 
of the algebra and they are given by aj = aj(J), ay G A. Such J's can be used to introduce 
ordering in the root system by assuming that a > if a( J) > 0. In what follows we will assume 
that all roots for which a(J) > are positive. Obviously one can consider the eigensubspaces 
of ad j as grading of the algebra g. 

In the case of symmetric spaces, the corresponding Cartan involution [33] provides a grading 
in g: = 0o © 0i> where 0o is the subalgebra of all elements of commuting with J. It contains 
the Cartan subalgebra t) and the subalgebra of g\t) spanned on those root subspaces g a , such 
that a(J) = 0. The set of all such roots is denoted by Ao- The corresponding symmetric space 
is spanned by all root subspaces in 0\0o- Note that one can always use a gauge transformation 
commuting with J to remove all components of the potential Q(x,t) that belong to 0o- 

For symmetric spaces of BD.I type, the potential has the form: 

f \ 

Q=\p s q , J = diag(l,0,...,0,-1). 
OjFso / 

For n = 2r + 1 the n-component vectors q and p have the form q = (qi, . . . , q r , qo, q f , . . . , qi) T , 



P = • • • ,Pr,Po,Pf > • • • >Pi) T ; while the matrix so = Sq enters in the definition of so 
X G so(n), X + S :n) X T S :n) = 0, where 



n): 



n+l 

s ( n) = E(- 1 ) s+1 < ) + i- 

s=l 

For n = 2r the n-component vectors q and p have the form q = (gi, . . . , q r , q F , . . . , qi) T , p 
(p 1 ,...,p r ,Pr,...,Pl) T and 

r 

: (n) _ V"^/ iu+l/ P (n) . rp(n) 



— / A~ L ) v s,n+l— s ~"~ n+l— s,s/ ' 



s=l 
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With this definition of orthogonality the Cartan subalgebra generators are represented by diago- 
nal matrices. By E^J above we mean nx n matrix whose matrix elements are (E^p)ij = 5 S iS p j. 

Let us comment briefly on the algebraic structure of the Lax pair, which is related to the 
symmetric space SO(n + 2)/(SO(n) x SO (2)). The element of the Cartan subalgebra J, which 
is dual to ei G E r allows us to introduce a grading in it: = 0o © 0i which satisfies: 

[x 1 ,x 2 ] g 00, [Xi.yi] g 01, [Yx,y 2 ] e 0O , 

for any choice of the elements Xi, X 2 G qq and Yi, Y 2 G 0i. The grading splits the set of positive 
roots of so(n) into two subsets A + = A^ U where A^ contains all the positive roots of 
which are orthogonal to e±, i.e. (a, ei) = 0; the roots in ft G A^ satisfy (ft, ex) = 1. For more 
details see the appendix below and [33] . 

The Lax pair can be considered in any representation of so(n), then the potential Q will take 
the form: 



Q( x ,t)= ^2 (q a (x,t)E a +p a (x,t)E. 



Next we introduce ra-component 'vectors' formed by the Weyl generators of so(n + 2) corre- 
sponding to the roots in Af: 

E l = (-E±(ei-e 2 )> ■ ■ ■ > E±( ei _ er ),E± ei ,E±( ei+er ^, . . . ,-E±( ei + e2 ))) 

for n = 2r + 1 and 

E t = ( E ±{e 1 -e 2 )i- ■ ■ !-S±(ei-e r );- E ±(ei+e r ) ) • • • > ^±(ei+e 2 ))> 

for n = 2r. Then the generic form of the potentials Q(x,t) related to these type of symmetric 
spaces can be written as sum of two "scalar" products 

Q{x,t) = (q(x,t) ■ Ef) + (p(x,t) • ). 

In terms of these notations the generic MNLS type equations connected to BD.I acquire the 
form 

nt + CLxx + 2 (Q,P)Q- (q,s q)s p = 0, 

ipt ~ Pxx ~ 2(q,p)p + (p, s p)s q = 0. (2.3) 

The Hamiltonian for the MNLS equations (|2.3h with the canonical reduction p = eg*, e = ±1 
imposed, is given by: 



#MNLS = J dx ({d x q, d x q*) - e(q, q*f + |(<f, s Q q)(q* , sq^ 5 )) 



2.1 Direct scattering problem for L 

The starting point for solving the direct and the inverse scattering problem (ISP) for L are the 
so-called Jost solutions, which are defined by their asymptotics (see, e.g. |18j and the references 
therein): 

lim (j)(x,t,\)e iXJx = 1, lim ^{x,t, X)e iXJx = t (2.4) 

X— > — oo x— too 

and the scattering matrix T(A) is defined by T(A,i) = ^^(^(x^t^X). Here we assume that the 
potential q(x,t) is tending to zero fast enough, when |x| — > oo. The special choice of J results 
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in the fact that the Jost solutions and the scattering matrix take values in the corresponding 
orthogonal Lie group SO(n + 2). One can use the following block-matrix structure of T(X,t) 



T(X,t) 



( ' m\ —b 7 

b + T 22 
\ c+ B+ T .< 



s B- 



T(X,t) 



I 



\ 



m 



i 

B~ 



B- T 

T22 

-b+ T 



sob' 



so rnf 



where ^(A, t) and B^(X, t) are n-component vectors, T22(A) and T22(A) are nxn block matrices, 
and m^(A), cf(X) are scalar functions. Here and below by 'hat' we will denote taking the 
inverse, i.e. T(X,t) = T _1 (A). Such parametrization is compatible with the generalized Gauss 
decompositions [33] of T(A,i). 

With this notations we introduce the generalized Gauss factors of T(A) as follows: 



T(X,t) 



TJ = e(^~ 



S+ = e (- + ^ + ) 



m 



_",± 



1 
mi 




^ s f- 
1 






D- 



1/mJ 



SJ = e(" 




=1 = S ° T )' 





->- T 


(• 


t 







1 





— T 


1 


it 


— T~' 



°1 
1 



s 1, 



where 



771 

m 



B+ 



' s p 



-:T 



m 



sot 



m, 



771V 



771 V 



mJ = T 22 + 



m. 



771 1 



T 22 + 



•SO 



b~b +T 



so 



771-, 



(2.5) 



These notations satisfy a number of relations which ensure that both T{X) and its inverse T{X) 
belong to the corresponding orthogonal group SO{n + 2) and that T(A)T(A) = IL. Some of them 
take the form: 



771+777^ + Q>-,B + ) + c£c[ = 1, 6+ B~ T + T22S0T22S0 + s B~b +J 's = 11 



1mt[cY 



b-' T sob~ 



0, 



77? 



T 22 B + - s B-cf = 0, 



2m i c\ 
mfB~ - 



B + > T s B + = 0, 



0. 



Important tools for reducing the ISP to a Riemann-Hilbert problem (RHP) are the fundamen- 
tal analytic solution (FAS) ^(x, i, A). Their construction is based on the generalized Gauss 
decomposition of T(A,i), see [HUTU [16]: 



X ± (x, t, X) = <P(x, t, X)Sf(t, A) = V(x, t, X)Tj(t, X)Dj(X). 



(2.6) 



More precisely, this construction ensures that £^(x,X) = x^{ x ^) el ^ Jx are analytic functions 
of A for A G C±. Here Sj, Tj 1 upper- and lower-block-triangular matrices, while Dj(X) are 
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block-diagonal matrices with the same block structure as T(X, t) above. Skipping the details, we 
give here the explicit expressions of the Gauss factors in terms of the matrix elements of T(A, t) 

S±(t, A) = exp ( ± (f±(A,t) • Ef)), T±(i, A) = exp ( =f (p*(A,t) • Ef)), 

where 

r f (A,t) = -^ T , f-(A,t) = — , p+(A,t) = -^, / T(A,t) = ^-, 
mj nil mJ 



and 



T 22 = mj - + , T 22 = m 2 — — — . 

T 22 = mj — r , T 22 = rn 2 



2m{ 2m l 

The two analyticity regions C+ and C_ are separated by the real line. The continuous spectrum 
of L fills in the real line and has multiplicity 2, see Section [3.3I below. 

If Q{x,t) evolves according to Q2.3j) then the scattering matrix and its elements satisfy the 
following linear evolution equations 

db^ t . o - . 1N dmf dm£ 

i— ±X 2 b ± (t,X) = 0, i^— ±X 2 B ± (t,X) = 0, i— ±- = 0, i—^- = 0, 
dt v ' 1 dt (it rfi 

so the block-diagonal matrices Z? ± (A) can be considered as generating functionals of the integrals 
of motion. It is well known |10| HI fT2j [14] that generic nonlinear evolution equations related 
to a simple Lie algebra g of rank r possess r series of integrals of motion in involution; thus 
for them one can prove complete integrability [3]. In our case we can consider as generating 
functionals of integrals of motion all (2r — l) 2 matrix elements of m^(A), as well as m^(A) for 
A £ C±. However they can not be all in involution. Such situation is characteristic for the 
superintegrable models. It is due to the degeneracy of the dispersion law of (j2.3j) . We remind 
that Dj(X) allow analytic extension for A G C± and that their zeroes and poles determine the 
discrete eigenvalues of L. 



2.2 Riemann— Hilbert problem and minimal set of scattering data for L 

The FAS for real A are linearly related 

X + (x,t,X)= X -(x,t,X)G ,j(X,t), G ,j(X,t) = Sj(X,t)Sj(X,t). (2.7) 

One can rewrite equation (|2.7p in an equivalent form for the FAS £^(x, t, A) = x^{ x i ^ X)e tXJx 
which satisfy the equation: 

i^ + Q{x)i ± {x,X)-X[J,£ ± {x 1 X)]=S 1 lim ^(x,t,X) = 1. (2.8) 

dx A->oo 

Then these FAS satisfy 

£+(x, t, A) = r (x, t, X)Gj(x, X, t), Gj(x, X, t) = e- iXJx G-j(X, t)e ajx . (2.9) 

Obviously the sewing function Gj(x, A, t) is uniquely determined by the Gauss factors Sj(X, t). 
Equation (|2.9|) is a Riemann-Hilbert problem (RHP) in multiplicative form. Since the Lax 
operator has no discrete eigenvalues, det^ ± (a;, A) have no zeroes for A G C± and ^(x,A) arc 
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regular solutions of the RHP. As it is well known the regular solution £ (x, A) of the RHP is 
uniquely determined. 

Given the solutions £ (x,t, A) one recovers Q(x,t) via the formula 

Q(x,t) = lim A(j-e ± Je ± (^,i,A)) = [J, £1(0;)], 

A— >oo 

which is obtained from equation (|2.8|) taking the limit A — >• 00. By £i(x) above we have denoted 
£i(x)= lim A(£(x,A)-l). 

A— s>oo 

If the potential Q(x, t) is such that the Lax operator L has no discrete eigenvalues, then the 
minimal set of scattering data is given by one of the sets Tj, i = 1,2 

Ti = {p+(A,t),p-(A,t), a G A+ AG R}, 

T 2 = {r+(A,t),r-(A ) t) ) a G A+ A G R}. (2.10) 

Any of these sets determines uniquely the scattering matrix T(A, t) and the corresponding po- 
tential Q(x,t). For more details, we refer to [THJ, [TT1 and references therein. 

Most of the known examples of MNLS on symmetric spaces are obtained after imposing the 
reduction: 

Q(x, t) = K^Q\x, t)K , K = exp ^ J^(3 + e 3 )Hj 

p k = K ql, K = diag (if ,22, • • • , #0,n+l,n+l) 

or in components p^ = eie^qt. As a consequence the corresponding MNLS takes the form: 

iqt + Qxx + 2(g, K q*)q- (q, s q)s K q* = 0. 

The scattering data are restricted by p~(X,t) = Kop*~'*(X,t) and r _ (A, t) = KoT + '*(X,t). 

If all ej = 1 then the reduction becomes the "canonical" one: Q(x, t) = Q'(x, t) and p~(X, t) = 
p^'*{X,t) and f"(A,t) =f+-*(A,t). 

2.3 Dressing factors and soliton solutions 

The main goal of the dressing method [S3 E31 ESI (Ml E21 CES] is, starting from a known solutions 
Xo(x,t, A) of Lo(A) with potential Q^{x,t) to construct new singular solutions x^(x,t, A) of L 
with a potential Q(i)(x,t) with two additional singularities located at prescribed positions A^; 
the reduction p = q* ensures that A^~ = (A^~)*. It is related to the regular one by a dressing 
factor u(x, t, A) 

xf (x,t, A) = u(x, A)xo (x,t, X)uZ 1 (X), u_(A) = lim u(x,X). (2.11) 

Note that u_(A) is a block-diagonal matrix. The dressing factor u(x, A) must satisfy the equation 
id x u + Q(i)(x)u - uQ(q){x) - A [J, u(x, A)] = 0, (2.12) 

and the normalization condition lim u(x,X) = 1. The construction of it(x, A) G SO(n + 2) is 

A— >oo 

based on an appropriate anzatz specifying explicitly the form of its A-dependence (see [48 [ [32] 
and the references therein). Here we will consider a special choice of dressing factors: 

u(x, A) = 11 + (c(A) - l)P(x) + f-L _ ^ P( x ), 
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P = S^P T S , c (A) = ^-|l (2.13) 

where P(x) and P(x) are mutually orthogonal projectors with rank 1. More specifically we have 

p(t f) - \n(x,t)){m(x,t)\ 

(n(x,t)\m(x,t)) ' { ' 

where (m(x,t)\ = (mo\(xo {x, X^)) -1 an< ^ \ n ( x ^)) = Xo~( x ' ^i~)l n o); ( m o| and |"o) are (constant) 
polarization vectors [23]. Taking the limit A — > oo in equation (|2.12p we get that 

Q (1) (x,t) - Q m {x,t) = (AJ- - X+)[J,P(x,t)-P(x,t)}. 

Below we assume that Qm\ = and impose Z 2 reduction condition 

KQK" 1 = Q, K = diag(ei,e 2 , . . . ,e r -i, l,e r ,e r _i, . . . ,e 2 ,ei). 

This in its turn leads to A^ = fi ± iv and \m a ) = K a \n a )*, a = 1, ... , n + 2. The polarization 
vectors (m(x, t)\ and |n(x,i)) are parameterised as follows: 

(m(x, t)\ = {m\{x, t),rri2(x, t), . . . , m r (x, t), 0, m r (x, t), . . . , m 2 (x, t),mi(x, t)) 

and 

|n(x,i)) = (ni(x,t),n 2 (x,t), . . . , n r (x, t), 0, n r (x, t), . . . , n 2 (x, t), n\(x, t)) T . 

As a result one gets: 

q[ ls) (x,t) = -2iv(P lk (x,t) + (-l) k P hn+2 (x,t)), 

where k = n + 3 — k. For more details on the soliton solutions satisfying the standard reduction 
Q(x,t) = Qt(x,t) see [291 [351 EH [211 [321 [22]. 

The effect of the dressing on the scattering data (|2.5p is as follows: 

^ b+ 1 ^ B~ 

Pi = — = -TTTPo > Pi = — = c(A)p ; 



mi 



Applying N times the dressing method, one gets a Lax operator with iV pairs of prescribed 
discrete eigenvalues A^, j = 1, . . . , N. The minimal sets of scattering data for the "dressed" Lax 
operator contains in addition the discrete eigenvalues A • , j = 1, . . . , N and the corresponding 
reflection/transmission coefficient at these points: 

Ti = {p^(A,t),p-(A,t),p+ )J (t),p+ j (t),A J ± , ae A+ Ael, j = l,...,N}, 

% = {r^(\t),T-(X,t),T^{t),T^(t),Xf, «GA+ A£t, j = l,...,iV}. (2.15) 



3 Resolvent and spectral decompositions 
in the typical representation of g ~ £? r 

Here we first formulate the interrelation between the 'naked' and dressed FAS. We will use 
these relations to determine the order of pole singularities of the resolvent which of course will 
influence the contribution of the discrete spectrum to the completeness relation. 

We first start with the simplest and more general case of the generalized Zakharov-Shabat 
system in which all eigenvalues of J are different. Then we discuss the additional construction 
necessary to treat cases when J has vanishing and/or equal eigenvalues. 
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3.1 The effect of dressing on the scattering data 

Let us first determine the effect of dressing on the Jost solutions with the simplest dressing 
factor ui(x, A). In what follows below we denote the Jost solutions corresponding to the regular 
solutions of RHP by ipo(x, A) and the dressed one by ip\(x, A). In order to preserve the definition 
in equation (|2.4[) we put: 



ipi(x, A) = m(x, X)ipo(x, A)«i j+ (A), 0i (x, A) = u\(x, \)4>o(x, A)«i _(A), 
where Ui±(x, A) = lim u\(x,X). We will also use the fact that iti±(A) are x-independent 

x— >±oo 

elements belonging to the Cartan subgroup of q. For the typical representation of so{2r + 1) 
and for the case in which only two singularities A^ are added we have: 



Ui +(A) = 1 + (ci(A) - l)E n + ( — !— - 1 ) E n+2 ,n+2, 



u(x, A) = exp (lnci(A)(Pi - Pi)) 

,ci(A) 

ui,_(A) = 1 + (ci(A) - l)K +2 ,n +2 + (^y - lj #11, 

«i ± (A) = exp(±lnci(A)J) . (3.1) 



Then from equation (|2.6p we get: 

xf(x,A) = ux(x, A)xo (x, A)«i _(A). 
As a consequence we find that 

Ti(A) = ui,+(A)T (A)«i _(A), Df(X) = u 1i+ (A)Aj(AK_(A), 

5f(A) = «i,_(A)5 (A)tli,_(A), r±(A) = «i j+ (A)r (A)ui,+ (A). (3.2) 

One can repeat the dressing procedure N times by using the dressing factor: 

u(x, A) = un(x, \)un-i(x, A) • • • ui(x, A). (3-3) 

Note that the projector P k of the fe-th dressing factor has the form of (|2.14p but the x-dependence 
of the polarization vectors is determined by the k — 1 dressed FAS: 

(x, A) = u k (x, A)ujfe_i(x, A) • • • ui(x, A)xo (x, A)ui _(A) • • • _(A)« fc _(A), 
u fc (x, A) = 11 + (c fe (A) - l)P fe (x) + f-i- - P fc (x), P fc = S^if S , (3.4) 

Cfe (A) = A ~ A fc t) = \Mx,t))(m k (x,t)\ 

Afc-A^' (n fe (x,f)|m fe (x,i)) ' 

(m k (x,t)\ = (m ,jfe|(Xfe_i(^,A^)) _1 , K(x,i)) = X j J"_ 1 (x, A^)|n ,fe), (3.5) 

and (mo,fc| and |roo,fc) are (constant) polarization vectors 

Using equation fl3.4f) and assuming that all Aj € C± are different we can treat the general 
case of RHP with 2N singular points. The corresponding relations between the 'naked' and 
dressed FAS are: 

T(A) = «_(A)r (A)u_(A), D ± (X) = u+(\)D (X)u-(\), 

S ± (X) = u-(X)S (\)u-(\), T±(A) = u + (X)T (X)u + (X), (3.6) 
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where 



U_(A) = t + (c(A) - l)£ n+2 ,n+2 + ( -^T - 1 ) #1,1, 



u±(A) = exp (± In c(A)J) , c(A) = J] c fc (A). (3.7) 

fc=i 

In components equations (|3.6p give: 

2 - m ro(A) 

m +(A) = mj (\)c (A), 174 (A) = ' , 
PtW = ^> P r(A) = c(AK(A), 

r! + (A) = ^, rf(A) = c(A)V(A), 

and m^(A) = mf Q (X). 

In what follows we will need the residues of u(x, \)xq(x, A) and its inverse Xoi x i tyu(x, A) 
at A = A^ respectively. From equations (I3.4p and (I3.5P we get: 

u(x, A)x+(s, A) * ^ T ~ 5y )(X) + X + ' (fe) W + 0(A - A+), 
u(s, A)xo (x, A) ~ (Afc+ " x Afc "^r W(x) + x~' (fe) (*) + 0(A - A^), 



(Afc 




) X +.(*)( X ) 




A- 


A^ 


(A fc + 


-A fc : 


) X -.(*)( X ) 




A- 


Afc 


(Afe 




)*+.(*) (x) 




A- 


A^ 


(A fc + 


-a,t; 


)£-.(*) (s) 




A- 


A* 



X +(x, A)u(x, A) ~ v fe - fc v y + A » + 0(A - A+), 

Xo(x, X)u(x, A) ~ v ' fe '- fc/ ;„ ^ + i" ik \x) + 0(A - Afc ), (3.8) 



where 





-,«(x 


x~ 




x H 


->«(x 


x~ 


->«(x 



X( fc _i)(a5, A fc )«fc+i(a;, A fc ) • • • tW(a:, A fc )P fc . (3.9) 
These results will be used below to find the residues of the resolvent at A = A^ • 

3.2 Spectral decompositions for the generalized Zakharov— Shabat system: 
s/(n)-case 

The FAS are the basic tool in constructing the spectral theory of the corresponding Lax operator. 
For the generic Lax operators related to the sl(n) algebras: 

Lgcn = i ~\~ (Qgcn(^) A t /g Cn )Xg Cn (x, A) — 0, {,Qgcn)jj{x) — 0j 
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this theory is well developed, see 07J 021 EJ- I- the generic case all eigenvalues of J gen = 
diag (Ji, J2, . . . , Jn+2) are different and non- vanishing: 

Jl > J2 > • • • > Jk > > Jfc+1 > • • • > J n+2 , tr J gcn = 0. 

The Jost solutions T/> gen (x, A), cj) gcn (x,X), the scattering matrix T gen (A) and the FAS Xgen^'A) 
are introduced by 03 HZ] (see also 021 El EH OB Q3B |22| ) : 

lim <f> gcn (x, \) e iJ ^ xX = 1, lim ip gcn (x, \) e u ^ xX = 1, 

X— > — 00 x— >oo 
^gen(A) = ^gen(x, A)0 gcn (x, A), 

XgenO^ A) = 0gcn(x, A)S± n (A), XgcnO^ A) = Vgen(x, A)T^ n (A)-Dj en (A), 

where (A), 7^ n (A) and D gcn {X) are the factors in the Gauss decompositions of T gon (A): 
Tgen(A) = r gen (A)D gen (A)5' geil (A) = r gen (A)D geil (A)5 gen (A). 

More specifically S^ cn (X) and T+ n (A) (resp. S gcn (X) and Tg^n(A)) are upper (resp. lower) trian- 
gular matrices whose diagonal elements are equal to 1. The diagonal matrices D gen (X) and 
D~ en (X) allow analytic extension in the upper and lower half planes respectively. 
The dressing factors u^x, A) are the simplest possible ones [37| : 

(-), 



U k , g en(x, A) = 1 + (cjfc(A) - l)Pfc(x), «fc )gen (~, A) = 1 + f - 1^ P k 



where the rank-1 projectors Pk{x) are expressed through the regular solutions analogously to 
equations (j3.3p ~ p.5p with Xo( x ,X k ) replaced by Xo, g en(-> K) 

The relations between the dressed and 'naked' scattering data are the same like in equa- 
tion ()3.2p only now the asymptotic values w ge n;±(A) are different. Assuming that all projectors Pk 
have rank 1 we get: 

A 1 



%on;±(A) = Ufe,gen;±(A), -fc,gen;±(A) = 1 + (cfc(A) - l)Pfc,±, 



k=l 

Pfc,+ = Es k ,s k , Pk,- = Ep k>Pk , 



where st (resp. pk) labels the position of the first (resp. the last) non-vanishing component of 
the polarization vector |~ofc). 

Using the FAS we introduce the resolvent P gen (A) of L gen in the form: 



R g en(X)f(x)= R gen (x,y,X)f(y). 

J —oo 

The kernel R gen (x, y, A) of the resolvent is given by: 



- (r x) _ / R ten& V, A) for A e C+, 

L genv ' ' 

where 



A) for A G 



^gen(-) 2/ J A) = ±ix gen (x, A)6 ± (x - y)x%n{Vi A), 

e ± (z) = (9(Tz)no-(9(± z )(i-no), n = ^p ss , 



s=l 
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Theorem 1. Let Q(x) be a potential of L which falls off fast enough for x — > ±00 and the 
corresponding RHP has a finite number of simple singularities at the points A^ G C±, i.e. 
X%n{ x i^) have simple poles and zeroes at . Then 

1) Rg Cn (x, y, A) is an analytic function of A for A G C± having pole singularities at G C±; 

2) R^ en (x, y, A) is a kernel of a bounded integral operator for Im A 7^ 0; 

3) R gen (x,y, X) is an uniformly bounded function for A G ffi and provides the kernel of an 
unbounded integral operator; 

4) R^ cn (x, y, A) satisfy the equation: 

L gcn (X)R± cn (x, y, A) = 18(x - y). 

Skipping the details (see [T7]) we will formulate below the completeness relation for the 
eigenfunctions of the Lax operator L gen . It is derived by applying the contour integration 
method (see e.g. [30l[l]) to the integral: 

d g cn(x,y) = -^—j> dXR gcn (x,y, A) - -^—^ j> dXR gcn (x, y, A), 

where the contours 7± are shown on the Fig. [1] and has the form: 

The explicit form of the dressing factors u gen (x, A) makes it obvious that the kernel of the 
resolvent has only simple poles at A = A^ . Therefore the final form of the completeness relation 
for the Jost solutions of L gen takes the form [T7] : 



n 

6(x-y)J2-E. 



ss 

1 a s 



— dX< 

2vr 



ko n 

J2\ X ^(x,X))(x^+(y,X)\- £ | X H-(*,A)><4t-(y,A)| 



+ ( Res R + (x,y) + Res R~(x,y) ) . 

\ A=A+ A=A~ / 

1=1 \ 1 J / 



(3.10) 
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It is easy to check that the residues in (|3.10p can be expressed by the properly normalized 
eigenfunctions of L gon corresponding to the eigenvalues A^ |17j . 

Thus we conclude that the continuous spectrum of L gen has multiplicity n and fills up the 
whole real axis R of the complex A-plane; the discrete eigenvalues of L gen constructed using the 
dressing factors u gen (x,X) are simple and the resolvent kernel R gen (x,y, A) has poles of order 
one at A = A^ . 

3.3 Resolvent and spectral decompositions for BD.I-type Lax operators 

In our case J has n vanishing eigenvalues which makes the problem more difficult. 
We can rewrite the Lax operator in the form: 

■°Xi . -sp -> x 

+ Q Xo = Axi, 

ox 

• <%) . . -, „ 

+ q xi + soqx-i = o, 

ox 

h<? 's Xo = Ax-i, 

ax 



where we have split the eigenfunction x(x, A) of L into three according to the natural block- 
matrix structure compatible with J: x(x,A) = (xi>Xo\x-l) ■ Note that the equation for xo 
can not be treated as eigenvalue equations; they can be formally integrated with: 

I'X 

Xo(x, A) = Xo,as + i dy (q*xi + s qX-l) , 

which eventually casts the Lax operator into the following integro-differential system with non- 
degenerate A dependence: 

i ~&x + iq ~ T ^ J dy (v* Xl + So ^~ 1 ) A ) = Axi, 

i ~^x~ + i ^^ s ° J dy (q* xi + So ^" 1 ) & A ) = -Ax-i- 

Similarly we can treat the operator which is adjoint to L whose FAS x( x i A) are the inverse to 
x(x, A), i.e. x(x, A) = x _1 ( x i A). Splitting each of the rows of x( x . A) into components as follows 
x(x,X) = (xi,Xo,X-i) we get: 

dxi 

i ~dx ~ ^°'^*' ) ~ Xxi = °' 

^ XiQ -X-i9'so = 0, 

ox 

■ dx 



- ~ (Xo,s q) - Ax-i = 0. 



dx 

Again the equation for xq can be formally integrated with: 

XoO, A) = xo,as + i j dy (xi(y, >^)f{y) + X-l(y, A)? t (y)s ) 



Now we get the following integro-differential system with non-degenerate A dependence 
i [ dy^xi{yA){f{y),q%x)) + x-i(yA){q\y)soq*(x))^+Xxi = 0, 



dx 
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■ dx 

i 



--i / dy ( K xi(y,^)(^(y)soq(x)) + x-i(yA)(q\y),q(x))) -Ax-i = o. 



dx 

Now we are ready to generalize the standard approach to the case of BD.I-type Lax operators. 
The kernel R(x, y, A) of the resolvent is given by: 



R(x,y,X) 

where 



J R + (x,y,X) for A G C+ 
\R (x,y,X) forAGC", 



R±(x, y, A) = ±i X ± (x, X)Q ± (x - ^(y, A), 

e ± (z) = e(Tz)En - e(±z)(t - e u ). 

The completeness relation for the eigenfunctions of the Lax operator L is derived by applying 
the contour integration method (see e.g. [3Q|. \T\) to the integral: 

3'(*,y) = ^-jf dAILi^x^A)-^- jf dAni^-^^A), 

where the contours 7± are shown on the Fig. Q] and II x = En +-E n +2,n+2- Using equations (|3.8|) 
and (|3.9p we are able to check that the kernel of the resolvent has poles of second order at 
A = A^ . Therefore the completeness relation takes the form: 

nif(*-y) = ^ /" dAn i {\x^(x,\))(x^(y,\)\ - \x [ ^ 2] -(x,X))(xit 2] ~(y,X)\} 

+ 2iY^ \ Res R + (x,y,X) + Res R~(x,y,X) > , 

where 

Res^^y, A) = ±(Afc - A+)nx (x + ' (fc) (x)^ + 'W(y) + x + ' (fe) (^)x + ' (fc) (y)) • 



A=A? 



In other words the continuous spectrum of L has multiplicity 2 and fills up the whole real axis M. 
of the complex A-plane; the discrete eigenvalues of L constructed using the dressing factors 
u(x, A) (|3.7p lead to second order poles of resolvent kernel R gen (x, y, A) at A = A^. 



4 Resolvent and spectral decompositions 
in the adjoint representation of g ~ B r 

The simplest realization of L in the adjoint representation is to make use of the adjoint action 
of Q(x) — XJ on q: 

£ a de a d = *~^T + IQ( X ) ~ ^ad, e^dix, A)] = 0. 

Note that the eigenfunctions of L ac j take values in the Lie algebra g. They are known also as the 
'squared solutions' of L and appear in a natural way in the analysis of the transform from the 
potential Q(x,t) to the scattering data of L, pQ; see also |14| 117] and the numerous references 
therein. 

The idea for the interpretation of the ISM as a generalized Fourier transform was launched 
in p]. It is based on the Wronskian relations which allow one to maps the potential Q(x, t) onto 
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the minimal sets of scattering data Tj . These ideas have been generalized also to the symmetric 
spaces, see [HI [2SJ [22J and the references therein. 

The 'squared solutions' that play the role of generalized exponentials are determined by the 
FAS and the Cartan-Weyl basis of the corresponding algebra as follows. First we introduce: 

e a,*i( x > A ) = X ± E a X ± (x, A), ef M (x, A) = X ± H j x ± (x, A), 

where x { X >X) are the FAS of L and E a , Hj form the Cartan-Weyl basis of q. Next we note 
that in the adjoint representation J ac [- = adj- = [J, •] has kernel. Just like in the previous section 
we have to project out that kernel, i.e. we need to introduce the projector: 

ttjX = adj 1 adjX, 

for any X £ q. In particular, choosing q ~ so(2r + 1) and J as in equation (|2.1j) we find that 
the potential Q provides a generic element of the image of 7rj, i.e. irjQ = Q. 

Next, the analysis of the Wronskian relations allows one to introduce two sets of squared 
solutions: 

=7rj( X ± (x,\)E aX ± (x,\)), ^t = ^j(x ± (x,X)E- aX ± (x,X)), aeA+. 

We remind that the set Ai contains all roots of so(r + 1) for which a(J) ^ 0. 
Let us introduce the sets of 'squared solutions': 

{*} = {*} c U {*} d , {*} = {*} c U{*} d , 

{*} c = {*+(x, A), *Z a (x, A), % < r, A e R} , 

f ■ + 1 N 

{*} d = *«.«,•(*), *- a;J W. *-«y( x )/. =1 ' 

{*} c = {*t a (x,A), *-(s,A), i<r, A6R}, 

{*} d = *^.(z), *« ;j (x)} =i , 



where the subscripts 'c' and 'd' refer to the continuous and discrete spectrum of L. 

Each of the above two sets are complete sets of functions in the space of allowed potentials. 
This fact can be proved by applying the contour integration method to the integral 



3g(x, y) = ^—j> d\G + (x, y, A) - ^- j> dXG (x, y, A), 

where the Green function is defined by: 

G ± (x, y, A) = Gf(x, y, X)6{y - x) - Gf(x, y, X)9(x - y), 
Gf(x,y,X)= £ A)® A), 

r 

Gf(x,y,X)= Yl ^± a {x,X)®^ a (y,X) + ^2hf(x,X)(g)hf(y,X), 

aGA UA~ i =1 

hf(x,X) = X ± (x,X)H jX ± (x,X). 
Skipping the details we give the result [25[ [27] : 
1 f°° 

5(x - y)U 0J = - / dX(G+(x, y, A) - G~(x, y, A)) 
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N 

-2iJ2(G+ j (x,y) + G^ j (x,y)), (4.1) 

3=1 

where 

n 0J = ( E <* ® ^-a - ® E <*), 

aGA+ 

4.1 Expansion over the 'squared solutions' 

The completeness relation of the 'squared solutions' allows one to expand any function over the 
'squared solutions' 

Using the Wronskian relations one can derive the expansions over the 'squared solutions' of 
two important functions. Skipping the calculational details we formulate the results [25] ■ The 
expansion of Q(x) over the systems {$ ± } and {^^j takes the form: 

• r-oo 

Q(x) = - / dA £ (r+(A)*+(x,A)-r-(A)*I a (x,A)) 

N 

+ 2 E E {<K,jW • • (4-2) 

■ /-oo 

Q(x) = — / dA ^ (p+(A)*l a (x,A)-p-(A)*-(x,A)) 

^ J —CO . i 

TV 

" 2 E E (pJj*^*) + ^*^(^)) • (4-3) 

The next expansion is of ady 1 5Q(.x) over the systems {^^j and {^ ,± }: 
i f°° 

BdfSQix) = — dA Y, (^+(A)*+(x,A) + «5r a "(A)*: Q (x,A)) 

J — OO i 

aeA+ 

JV 

+ E E ( 5W U x ) - 5 ' w -^)) > (4-4) 

i f 00 

ad7 1 5Q(x) = — ^ dA Y, W(A)*i Q (x,A) + ^-(A)*-(x,A)) 

~°° aGA+ 

JV 

+ E E { 6W -a-^) ~ SW^to) , (4-5) 

fc=l aeA+ 

where 

^± a;i (x) = a±r Q ±.*± Q;J (x) + 5r Q ±.*± a;J (x), 
and $% a;j (x) = &%Jx,\f), *± Qy (x) = a A *^ Q (2;, A)| A=A ±. 
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The expansions (|4.2p . (|4,3p is another way to establish the one-to-one correspondence bet- 
ween Q(x) and each of the minimal sets of scattering data 7\ and T2 (|2.1(jp . Likewise the 
expansions (|4.4p . (|4.5p establish the one-to-one correspondence between the variation of the 
potential 5Q(x) and the variations of the scattering data 57± and 872- 

The expansions (j4.4|) . (|4.5j) have a special particular case when one considers the class of 
variations of Q(x,t) due to the evolution in t. Then 

SQ(x, t) = Q(x, t + St)- Q{x, t) = ^5t + 0({5t) 2 ) . 

Assuming that 5t is small and keeping only the first order terms in 5t we get the expansions for 
adj 1 Qf. They are obtained from (|4.4p . (|4.5p by replacing 5p^(\) and 5t^(\) by dtp^{X) and 



4.2 The generating operators 

To complete the analogy between the standard Fourier transform and the expansions over the 
'squared solutions' we need the analogs of the operator Dq = — id/dx. The operator Dq is the 
one for which e lXx is an eigenfunction: Doe lXx = Xe lXx . Therefore it is natural to introduce the 
generating operators A± through: 



(A+ - X)*± a (x, A) = 0, (A+ - A)*"(x, A) = 0, 
(A_ - A)*+(x, A) = 0, (A_ - A)*I a (x, A) = 0, 

where the generating operators A± are given by: 



(A + -\f)*+ a . j (x) = 0, 
(A + -\f)*+ a . j (x) = 0, 



dX 



dx 



Q{x) 



dy[Q(y),X(y)] 



±00 



The rest of the squared solutions are not eigenfunctions of neither A + nor A_: 



(A + -\+)4>Z a . j (x) = *± a . 
(A_-A+)*i ;i (x) = *+ ;i (a 



», (A + — \j )^f a ,j(x) = * ay .(x) 
(A_ - Aj)*~ w -(x) = *~j-(x) 



i.e., and are adjoint eigenfunctions of A + and A_. This means that A - , j = 

1,...,N are also the discrete eigenvalues of A± but the corresponding eigenspaces of A± have 
double the dimensions of the ones of L; now they are spanned by both ^^ a .j(x) and f ^(s). 
Thus the sets {^/} and {$} are the complete sets of eigen- and adjoint functions of A + and A_. 

Therefore the completeness relation (|4.ip can be viewed as the spectral decompositions of the 
recursion operators A±. It is also obvious that the continuous spectrum of these operators fills 
up the real axis M. and has multiplicity 2n; the discrete spectrum consists of the eigenvalues A^ 
and each of them has multiplicity 2. 



5 Resolvent and spectral decompositions 
in the spinor representation of g ~ B r 

Using the general theory one can calculate the explicit form of the Cartan-Weyl basis in the 
spinor representations. In Appendices iBl and ICl below we give the results for r = 2 and r = 3. 
Therefore in the spinor representation the Lax operators take the form: 

dips 

L S pil) S p = i-^- + (Qs P ~ \J S p)ip sp (x, A) = 0, (5.1) 
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where Q sp (x,t) and J sp are 2 r x 2 r matrices of the form: 

^ sp ~ W o) ' sp " 2 V o -t 2 ) ' 

where the explicit form of q{x) for r = 2 and r = 3 is given by equations (jB.ip and (jC.ip below. 

It is well known [5] that the spinor representations of so{2r + 1) are realized by symplectic 
(resp. orthogonal) matrices if r(r + l)/2 is odd (resp. even). Thus in what follows we will 
view the spinor representations of so(2r + 1) as typical representations of sp(2 r ) (resp. so(2 r )) 
algebra. Combined with the corresponding value of J sp (|B.1|) one can conclude that that the 
Lax operator L sp for odd values of r (r + l)/2 can be related to the C. Ill-type symmetric spaces. 
The potential Q sp (x) however is not a generic one; it may be obtained from the generic potential 
as a special reduction, which picks up so(2r + 1) as the subalgebra of sp(2 r ). Below we will 
construct an automorphism whose kernel will pick up so(2r + 1) as a subalgebra of sp(2 r ). 

Similarly the Lax operator L sp above for even values of r(r + l)/2 can be related to the 
so(2 r ) algebra. The element J sp (|B.1|) is characteristic for the D. Ill-type symmetric spaces. 
The potential Q sp (x) may be obtained from the generic potential as a special reduction, which 
picks up so(2r + 1) as the subalgebra of so(2 r ). 

The spectral problem (|5.ip is technically more simple to treat. It has the form of block-matrix 
AKNS which means that the corresponding Jost solutions and FAS are determined as follows: 

if>(x,X) ~ e~ lXJx , cf>(x,X) ^ e~ iXJ *, T(X) = ( £ ) , 

X— >oo X^f — OO \ O CL J 

iP(x,\) = (ip~(x, X),tp + (x,X)), 0(x,A) = (0 + (x,A),</»-(x,A)), 
X + (x,A) = (0+(x,A),V + (x,A)), X '(x,X) = (^-(x,A),0-(x,A)). 

5.1 The Gauss factors in the spinor representation 

The spectral theory of the Lax operators related to the symmetric spaces of C.III and D.III 
types were developed in |2CH [28] . What is different here is the special choice of the rank and 
the additional reduction 7TB r which picks up the spinor representation of so(2r + 1). In our 
considerations below we will assume that this reduction is applied. So though all our 2 r x 2 r 
matrices are split into blocks of dimension 2 r_1 x 2 r ~ 1 , the corresponding group (resp. algebraic) 
elements belong to the (spinor representation of) group SO(2r + 1) (resp. algebra so(2r + 1)). 
Thus we define the FAS of L sp by: 

X t p (x,X) = (|<^ + ), |^ + c+)) (x, A) = <f>(x,X)S+ (A) = A p (x, A)T-(A)£+(A), 

A) = (|^-c->, Ic/O) (x,A) = <t>(x,X)S- p (X) = i> sp (x,X)T+(X)D-(X), (5.2) 

where the block-triangular functions S^A) and T^ p (X) are given by: 
S s + p(A) = ( o ^ \ (A) ) , T~ p (A) = ( h+ l+ {X) I ) , 

wi^- w :)' riw-(;-*T w ). (5.3) 

The matrices D^p(X) are block-diagonal and equal: 

^) = ( O+ (A) 4)). O i (*)-(T ) .- W )- 

The supper scripts ± here refer to their analyticity properties for A E C±. 
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All factors S^ p , T^ p and D^p take values in the spinor representation of the group SO(2r + 1) 
and are determined by the minimal sets of scattering data (|2.15|) . Besides, since 

T sp (A) = T- p (\)D+(\)st p (\) = T+(X)D-(X)S; p (X), 
f sp (A) = S+ (A)£+ (A)T~ (A) = S-(\)D-{\)T+{\), 

we can view the factors Sf p , T^ p and D^p as generalized Gauss decompositions (see [53]) of T sp (A) 
and its inverse. 

From equations (|5.2|) . (|5.3p one can derive: 

xt P ( x , A) = X7 P (x, A)G , sp (A), Xs P ( x > A) = xtp(x, A)Go, sp (A) ; (5.4) 

G o,sp(A) = 1 + r r + _ r+ ) , G ,sp(A) = ( 1 + _ T t + J~ ~l + ) (5.5) 

valid for A G R. Below we introduce: 

X±(x,\)=xt P ( x ^> iXJX - ( 5 - 6 ) 



'±1 

sp> 

nice property, namely their asymptotic behavior for A — > ±oo is given by: 



Strictly speaking it is X^(x, A) that allow analytic extension for A G C±. They have also another 



lim X±(x,\) = H.. (5.7) 

A— >oo 

Along with X^(x, A) we can use another set of FAS X^(x, A) = X^(x, \)D^p, which also satisfy 
equation (|5.7|) due to the fact that: 

hm Df p (X) = 1. 

The equations (|5.4p and (|5.5p can be written down as: 

X+(x, A) = X sp (x, A)G sp (x, A), A G R (5.8) 

with 



G sp (z, A) = e- iA ^G 0! s P (A)e iA ^, G , sp (A) = ( 



11 + p~p + p~ 



Equations (|5.8p combined with (|5.7p are known as a Riemann-Hilbert problem (RHP) with ca- 
nonical normalization |13j . It has unique regular solution; the matrix- valued solutions Xq (x, A) 
and Xq sp (x, A) of (|5.8p . (|5.7p is called regular if det X^ sp (x, A) does not vanish for any A G C±. 
One can derive the following integral decomposition for X^(x, A): 

v 



X S >,A) = 1 + ^ /I ^X~(x, p) Kl (x, p) + Yl _ 1 ■ M 

1 f°° da Xf(x)K 2 ,j{x) 
X-(x,\) = II + — J ^ -^- X -(x,»)K 2 (x,ri - £ A+ _ A ( 5 - 10 ) 



where Xf (x) = X^(x,\f) and 
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Equations (|5.9p , (|5.10p can be viewed as a set of singular integral equations which are equivalent 
to the RHP. For the MNLS these were first derived in [3D] . 

Finally, the potential Q(x,t) can be recovered from the solutions X^(x, A) of RHP (|5.8p with 
a canonical normalisation (|5.7p . Skipping the details, we provide here only the final result: 

Q sp (x,t) = lim \(J-X±(x,\)JX±(x,\)]) = [J,Xx(x)], 

A->oo 

where X\(x) = lim (X sp (x, A) — 1). 

A— >oo 

5.2 Definition and properties of R^(x, y, A) 

The resolvent P sp (A) of L sp is again expressed through the FAS in the form: 



/oo 
R sp (x,y,X)f(y), 
-oo 

where R sp (x,y,X) are given by: 

R+(x,y,X) forAGC+, 



s P v , l^-^^A) for AG 



and 



Rf p (x, y ,X) = ±i x f p (x,X)Q ± (x-y)xf P (yA), e ± (z)=( ^ )t _ e{ l z)1 

5.3 The dressing factors in the spinor representation 

The asymptotics of the dressing factor for x — > ±oo are: 

u± = exp (lnci(A)J S p) , 

where 



,^ 2 r \ _(l/^t 2 r-l 

+ V 1/^2-1 J ' V \/^ I1 2-l 

Thus the asymptotics of the projectors P and P in the spinor representation are projectors of 
rank 2 r_1 . Since the dynamics of the MNLS does not change the rank of the projectors we 
conclude that the dressing factor in the spinor representation must be of the form: 

u(x,X) = exp f ilnci(A)(Pi - Pi) ) = V / ^A)F t (M) + ?=t A(M) 

/ v c n A ) 

dressing factors of such form with non-rational dependence on A were considered for the first 
time in by Ivanov [35] for the MNLS related to symplectic algebras. Here we see that such 
construction can be used also for the orthogonal algebras. 

However when it comes to analyze the relations between the 'naked' and the dressed scattering 
matrices and their Gauss factors we get integer powers of c±(X). Indeed, for the simplest case 
when the dressing procedure is applied just once we have: 

T sp (A) = n + T 0iSp (A)^(A), D±(A) = u + D± sp (X)u„(X), 
T±(A) = n + T ± sp (A)u + (A), S±(A) = n_5 ± sp (A)n_(A), 
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and as a consequence 

a t P W = ^j^y a <ts P ( A )> a s~ P ( A ) = 4( A ) a ^ sp ( A ) 5 

ptpW = ^+^yP^ sp ( A ). Psp( A ) = c i"( A Hs P ( A )' 

T sp( A ) = ^X) T otsp( A )' T sp( A ) = 4( A ) T oTsp( A )- 

5.4 The spectral decompositions of L sp 

Again apply the contour integration method to the kernel R^p and derive the following com- 
pleteness relation: 

1 f°° 

8(x-y)Z 2 r = — J dA{|0 + (x,A))a + (A)(V + (y,A)| -|0-(x,A))a-(A)(V^(y,A)|} 

+ V files R + (x,y) + Res R~(x,y)\ . (5.11) 

^- \ \A=A + \=\7 

3=1 \ 3 J / 

The residues in (|5.1ip can be expressed by the properly normalized eigenfunctions of L sp corre- 
sponding to the eigenvalues A^ . 

Thus we conclude that the continuous spectrum of L sp fills up M and has multiplicity 2 r . 
The discrete eigenvalues A^ are simple poles of the resolvent. 

6 Dressing factors and higher representations of q 

Here we will briefly outline how, starting from equation (I2.13|) one can construct the dressing 
factors in any irreducible representation of the Lie algebra q. We will illustrate this on one of 
the simple nontrivial examples of u(x, A) with rank-1 projectors P{x) and P{x). Our intention 
is to outline the explicit A-dependence of u(x, A) in any IRREP. To this end it will be most 
convenient to use the first line of equation (|3,1|) : 

u(x, A) = exp (lnci(A)(Pi - Pi)) ■ (6.1) 

Note that P(x) — P(x) £ g and therefore the right hand side of (16. ip will be an a Lie group ele- 
ment. In what follows we will assume that rankP(x) = rankP(x) = 1 and will derive explicitly 
the A-dependence of the right hand side of equation (|6.1h in any irreducible representation of g. 
In fact it will be enough to analyze the A-dependence of the asymptotic of u(x, A) for x — > ±oo. 

lim (P(x) - P») = H ei . 

x— >oo 

In order to be more specific we will do our considerations for the case when f ~ so(7). This 
algebra is of rank 3. We also choose the representation with highest weight 

oj = 3w 3 = -(ei + e 2 + e 3 ). 

The structure of the weight system is described in Table [TJ 

It is natural to expect that u(x, A) will have the same type of A-dependence [15] as its 
asymptotic for x — > ±oo. Therefore we will evaluate the right hand side of equation (|6.ip for 
x — > oo. Doing this we well use the known formula 



H ei = ( e i<7)l7)<7l- 
7 er( 3 "3) 
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Table 1. The structure of the weight system r^ 1 ) for the algebra so(7) with dimension 112. We list 
here the number of weights of different lengths and their multiplicity ^(7) and length. The indices 
i, j, k are different and take values 1, 2 and 3. 



weight type 




r 2 


r 3 


r 4 




|(±ei±e 2 ±e 3 ) 


^(±3a ± 3ej ± ek) 


l(±3ei ± ej ± ejfc) 


|(±ei ± e 2 ± e 3 ) 


#(7) 


8 


24 


24 


24 




1 


1 


2 


4 


%) 


27 
4 


19 
4 


11 
4 


3 
4 



Therefore we have to arrange the weights in I^ 3 " 3 ) according to their scalar products with e±, 
namely: 

r (3tJ3) = r 3/2 u r 1/2 u r_ 1/2 u r_ 3/2 , 

T 3/2 = ||(e x ± e 2 ± e 3 )| U ji(3ei ± 3 ej ± e fc ) J U ji(3ei ± e 2 ± e 3 )| , 



1/2 



-(ei ± 3e 2 ± 3e 3 ) 1 U j -{e x ± 3 ej ± e fc ) 1 U | -(ei ± e 2 ± e 3 ) 



r-i/ 2 = |^(-ei ± 3e 2 ± 3e 3 )| U |^(-ei ± 3 ej ± e fc )| U |^(-ei ± e 2 ± e 3 )| 

r_ 3/2 = {|(-ei ± e 2 ± e 3 )| U j^(-3ei ± 3e^ ± e fc )| U j^(-3ei ± e 2 ± e 3 ] 

where j, k take the values 2 and 3. As a result we obtain 
3 1 1 3 

H ei = 2 7F 3/2 + 2^1/2 ~ 2 7r ~ 1 /2 ~ o 71 "- 3 / 2 ' 

where the projectors 7r a , a = ±§, ±~ are equal to 



(6.2) 



and obviously satisfy the relations: 

TTaVTb = (5 ab 7r a , rank7r 3 / 2 = rank7r_ 3 / 2 = 20, 



rank7Ti/ 2 = rank7r_i/ 2 = 36, 



and 7r 3 / 2 + 7Ti/ 2 + vr_ 1 / 2 + 7r„ 3 / 2 = Un 2 . Thus the Q{x) — A J acquires the following block-matrix 
form: 

/ f Al 20 Q(12) Q(i3) Q(u) ^ 
Q(21) 2 A1 3 6 Q(23) Q(24) 
Q(31) Q(32) -^Al 36 Q( 34 ) 
V Q(41) Q(42) Q(43) — 1-^120/ 

Formally this potential can be viewed as related to the homogeneous space S'O(112)/5(O(40) ® 
0(72)). However all matrix elements of the potential Q(x) are determined by the five com- 
ponents of the vector q and their complex conjugate. This deep reduction imposed on Q(x) 
corresponds to the fact that instead of considering a generic element of this homogeneous space, 
we rather pick up the representation of so(7) with highest weight 3u; 3 . 
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Inserting equation (|6.2p into equation (|6.ip we get: 

lkn u^\x,X) = (c(X)f 2 7r 3/2 + WA)) 1 ^^ + ( c (A))^ 1 /2 7r _ i/2 + (^A))" 3 / 2 ^/,, 

and as a result for the A-dependence of tt^ 3w3 ^(x, A) and its inverse we get: 

u^\x,X) = (c(X)f/\ 3/2 (x) + (c(A)) 1 /2 vri/2 (x) 

+ (c(A))" 1 / 2 vr_ 1/2 ( a ;) + (c(A))- 3 / 2 7r„ 3/2 (x), 
(u^rHx, A) = (c(A))" 3 / 2 vr 3/2 ( 3 ;) + (c(A))- 1 / 2 vr 1/2 (x) 

+ ( C (A)) 1 / 2 vr„ 1/2 (x) + (c(A)) 3 / 2 vr_3/ 2 (x). 

The four projectors 7T a (x) have the same properties as their asymptotic values: 

TTa(x)ir b (x) = S ab TT a (x), 

rank -^3^(2;) = rank-7r_3/ 2 (x) = 20, rank-/r 1 / 2 (x) = rank7r_ 1/ / 2 (2;) = 36. 

Their explicit x dependence as well as the interrelation between the potentials Q(o)(x) and 
Qn)(x) follow from the equation for u(x, A) f|2. 12|) considered in the representation y( 3w3 ). In 
particular we get: 

Q{i){x) - Q(p)(x) = lim A (j - u^Jiu^y 1 ^, A)) 

3 1 

= (A + - A") J, -(7T 3/2 (x) - 7T„ 3/2 (x)) + -(7T 1/2 (x) - TV-i/ 2 (x)) 

Though the expressions for the dressing factor in this representation seem to be rather complex, 
nevertheless they are determined uniquely by the projector P\(x), or equivalently, through the 
polarization vector \ni(x)). The corresponding expressions can be using the fact that can 
be extracted as the invariant subspace of the tensor product V^ 3 ^ ® V^ 3 ^ <£> V^ 3 ^ corresponding 
to the highest weight vector 3^3. Therefore one can conclude that the matrix elements of the 
projectors 7r a , a = ±3/2, ±1/2 will be polynomials of sixth order of the components of |ni(x)). 

Our final remark concerns the analyticity properties of the FAS dressed by u^ Sul3 \x,X). 
The factor itself contains powers of root square of c(A) which in general may lead to essential 
singularities at A = A + and A = A~. Note however that the dressed FAS in this representation 
are obtained from the regular solutions via the analog of equation (|2. 1 1 [) : 

xf' 3 " 3 (x,t,X) = (x, X) X ^ 3 (x,t, A) (u 3 ^ 3 )- 1 ^) , 
u^ 3 (X) = lim n 3aJ3 (x,A) =exp(-lnci(A)F ei ) 

X— > — 00 

= (c(A)) 3 / 2 vr =3/2 + (c(A)) 1 / 2 7r_ 1 / 2 + {c{X)r l '^ 1/2 + ( C (A))" 3 / 2 vr 3/2 . (6.3) 



It is not difficult to check that the right hand side of first line in equation (I6.3P does not contain 
square root terms of c(A); all half-integer powers of c(A) get multiplied by other half- integer 
powers and the result is that Xi' 3uJ3 ( x i t, A) acquires additional pole singularities at A = A + and 
A = A~. 



7 Conclusions 

We have analyzed the spectral properties of the Lax operators related to three different repre- 
sentations of q ~ B r : the typical, the adjoint and the spinor representation. In all these cases 
the spectral properties such as: i) the multiplicity of the continuous spectra and of the discrete 
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eigenvalues; ii) the explicit form of the dressing factors; iii) the completeness relations of the 
eigenfunctions are substantially different. However the minimal sets of scattering data Tj are 
provided by the same sets of functions, i.e. the sets Xj are invariant with respect to the choice 
of the representation of Q. 

Our considerations were performed for the class of smooth potentials Q(x) vanishing fast 
enough for x — > ±00. Similar results can be derived also for the class of potentials tending to 
constants Q± for x —> ±00. 

A The adjoint representations of so(5) and so(7) 

The adjoint representations of all orthogonal algebras so(2r + 1) and so(2r) are characterized 
by the fundamental weight w 2 = &i + e 2 - By definition the corresponding weight system is 

r W2 = Au{o} r , 

where A is the root system of the algebra and {0} is a vanishing weight with multiplicity r. We 
will order the weights in r w2 according to their scalar products with e±. 

Let us consider in more detail the two special cases of so(2r + 1) with r = 2 and r = 3. For 
r = 2 the adjoint representation is 10-dimensional. Ordering the roots as mentioned above we 
get: 

&t ~ {ei + e 2 ,ei,ei - e 2 }, A ~{e 2 ,-e 2 }, Aj" ~ {-ei + e 2 , ei, -&\ - e 2 }, 
As a result the element J in the adjoint representation takes the form: 
/1 2 \ / Q ad; i 2 

Jad = I 06 J , Qad = I Qad;21 06 Qad;23 

\ -1 2 J V Qad;32 

Analogously for r = 3 the adjoint representation is 21-dimensional. Ordering the roots as 
mentioned above we get: 

A^ ~ {ei + e 2 ,ei + e 3 ,e 1( ei - e 3 ,ei - e 2 }, A ~ {e 2 ± e 3 , -(e 2 ± e 3 )}, 
A r - { _e i + e 2, -ei + e 3 , ei, -ei - e 3 , -ei - e 2 }. 

As a result the element J in the adjoint representation takes the form: 

/H 5 \ / Qad;12 

4d = On I , Qad = I Qad;21 On Q a d;23 

V -1 5 / V Qad;32 

It is not difficult to write down the explicit form of Q a d but it will not be necessary. We have 
used above a more compact realization of Q a d as adQ. 

B The spinor representation of so(5) 

The highest weight and the weight system of so(5) are given by [5j [33]. It is well known that 
so (5) — sp(4) so the spinor representation of so(5) is realized through symplectic sp{A) matrices 

!/ 1/ 
^2 = 71 = + e 2j, 72=2V e i _e 2)' 73 = -72 , 74 = -7i- 
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The Cartan-Weyl basis of so(5) is given by 

E ei -e 2 = r 2j 2 = £2e 2 j -^61+62 = ^1,1 = &2ei, 

Eei = ^1,2 ~ ^2,1 = £<Ei+e 2 ' ^e 2 = ^1,2 — T 2 J = £ ei _ £2 , 

1 1 

where k = 5 — k and 

Fk, P = |7fe)(7pl> l<fc<p<4. 

By £ ei ±ej above we have denoted the Weyl generators of sp(4). The Lax operator L in the spinor 
representation of so(5) takes the form: 

dips 

Lspfpsp = i-r^r + (Qsp - AJ sp )^ sp (a;, A) = 0, 
where Q sp (x,t) and J sp are 4x4 symplectic matrices of the form: 

«-=(?.?)• ^Ko'-O- ^=(S4)- (Bi) 



C The spinor representation of so(7) 

The highest weight and the weight system of so(7) are given by [5j [33] 
^3 = 7i = 2 ( ei + e 2 + e 3j, 72 = ^(ei + e 2 - e 3 ), 

73 = ~ e 2 + e 3j, 74 = ~ e 2 - e 3), 

75 = -74, 76 = -73 , 77 = -72 , 78 = -71- 

The Cartan-Weyl basis of so(7) is given by 



E ei —e 2 


— ^3,4 




E ei - e:j 


= ^2,4 


= ^^2+64) 




= r a 





Ee2 — e-i 


— ^2,3 


— p 

^£2 


-£3 ! 




= r 1>3 




f£4) 




= r M 


= e £1 


-64 5 



^ei — I\ 4 + r 2) 3 — £ ei +e 4 + £ £2 +e3> E e2 — Tl^ + 1^2,4 — £<Ei-e 4 + £e2-£4> 

E e3 = ^1,2 — T3,4 = £ei-e 2 — ^e 3 -e 4 , H ei = nC^Yl + ^2,2 + ^3 + ^4), 

= 2^ 1 ' 1 ^2,2 — ^3 — r^), H ei = 2(^1,1 ~~ ^2,2 + r3 ? 3 — r^), 

where k = 9 — k and 

r k , P = |7 fc )(7 P l - (-l) fc+p l7 P -)(7fcl, 1 < k < p < 4, 
r k ,p = |7fc)(7 P l + (-l) fc+p l7 P -)(7fcl, 1 < *: < P < 4. 

Note that the typical representation of so(8) is also 8-dimensional. So by £ £i ±e- above we have 
denoted the Weyl generators of so(8). So we can also consider the spinor representation of so(7) 
as an embedding of so{7) into so{8). Such embedding can be realized by taking the average of the 
Cartan-Weyl basis of so(8) with respect to the Z2 external automorphism of so(8) (the mirror 
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reflection) which changes €4 f-> —64. In short the Lax operator L in the spinor representation of 
so(7) take the form: 

dips 

Lsptpsp = i-g^r + (Qs P - XJ sp )ip sp (x, A) = 0, 
where Q sp (x,t) and J sp are 8x8 matrices of the form: 

o-p = (°tj). ^ = \{ o* -u ) ' q{x ' t] 

Our final remark here concerns the spinor representation with highest weight u = Suj^. Then 
f 3 1 

r ~ <^ o^ ±ei ± 62 ^ e3 ^' 2 ^ ±ei ± 62 ± e ^ 

The dimension of this representation is 8. The explicit form of the element J in this representa- 
tion is J = diag ^1, — |l, — i.e. all eigenvalues of J are non-vanishing. The potential Q 
will have block-matrix structure compatible with the one of J. 



( 90 95 9l \ 

92 9o -<?i 
\ -qi q 2 -90 / 



(C.l) 
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